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O ' Abstract 



We show that lepton number violating muon decays, /x + —* e + + + V\ (i = 
e,fi or r), can consistently explain the neutrino anomaly reported by the LSND 
experiment. Two effective operators in the Standard Model are identified which 
lead to just such decays and no other processes. The scale of new physics A must 
be relatively low, A < 500 GeV. Extensions of the Standard Model which realize 
these effective operators are presented. Since new physics affects only the decay of 



the muon, and not of 7r ± , our scenario predicts a null result for Un — v e oscillation 
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searches at the Fermilab mini-BOONE experiment. Models which realize these 
effective operators, while consistent with all available data, can be tested in the 
near future through (i) discovery of new scalar particles with masses below about 
500 GeV, (ii) small but observable deviations in e + e~ — > and e + e~ — > uv^y 

cross sections, (iii) observable corrections to the muon g— 2, and (iv) lepton number 
violating Z° decays with branching ratios of order 10 -7 . 



1 Introduction 



Evidence for neutrino flavor oscillations have been mounting over the years. A variety 
of solar neutrino experiments [I], @, [§, |3J, §] which have consistently detected fewer z/ e 's 
from the sun than expected now seem to converge on the large angle MSW solution as 
the preferred explanation for the discrepancy ||. The deficit in the flux of atmospheric 
muon neutrinos J7| and especially the observed zenith angle dependence are compelling 
evidences in favor of z/ M — v T oscillations. The LSND collaboration observes an anomaly 
in the flux of V e detected ||, which can be interpreted as evidence for V^—Ve oscillations. 

As it turns out, it is not possible to explain all of these observations in terms of neu- 
trino oscillations with just three flavors of neutrinos (V e , and z/ T ). The reason is that 
the characteristic oscillation length scales inferred from the three types of experiments 
are quite different having no overlap: A ~ 200 km, X atm ~ 6000 km and Xlsnd ~ 60 
m. The oscillation length is given by A = (AE)/Am 2 where Amf 2 = m\ — m\ etc, 
with rrti being the mass of the ith neutrino mass eigenstate, and where E denotes the 
neutrino energy. In three neutrino oscillation schemes there are only two independent 
mass-splittings, which therefore excludes a simultaneous explanation of all three obser- 
vations. 

Two approaches have been adopted in the literature to address this conundrum. In 
the first, one of the three observations is simply discarded. In the second approach, a 
fourth "sterile" neutrino (u s ) - sterile so that it does not affect the number of neutrinos 
that couples to the Z° bososn - is introduced. Both of these approaches are not very 
satisfactory. The first one is clearly without justification. The second, introduction of a 
sterile neutrino, has some serious theoretical difficulties, in addition to a possible prob- 
lem with the standard Big Bang Nucleosynthesis. Its lightness cannot be explained in 
any simple way, unlike that of the "active" neutrino which follows very naturally from 
the seesaw mechanism 0. Even taking a light sterile neutrino for granted, some po- 
tential experimental difficulties have emerged over the last year. When the recent SNO 
results on solar neutrinos are combined with the results of SuperKamiokande, the al- 
lowed parameter space is tightly constrained for oscillations involving a sterile neutrino. 
In fact, a two-flavor v e — v s oscillation scenario no longer provides a good fit to the 
solar neutrino data. Atmospheric neutrino oscillation results from SuperKamiokande 
also disfavor oscillations of with a v s . Somewhat more involved four neutrino oscilla- 
tion schemes would be necessary, conforming to either a (3+1) |lJJ] scheme or a mixed 
2+2 scheme [ITT]. Although both appear to be viable presently, neither gives a really 



satisfactory global fit to all the data [12]. There is also a proposal to account for the 



LSND results in the three neutrino framework by invoking CPT violation [ IB ]. 

The purpose of this paper is to provide a simultaneous solution to the three neutrino 
anomalies without introducing a sterile neutrino. It will be achieved by small non- 
standard interactions of the leptons instead. Since the flavor conversion probabilities 
are of order unity for both the solar and the atmospheric neutrino experiments, we shall 
focus on non-standard lepton interactions explaining the LSND anomaly, which calls for 
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a smaller probability: P v ^v e = (0.264 ± 0.067 ± 0.045)% [g. Specifically, we shall show 
that lepton number violating (AL = 2) muon decays can account for the LSND events. 
The new decay mode of the muon is 

H + -> e + + V e + Vi , (1) 

where Vi stands for any one of z/ e , z/ M or V T . Consistent gauge models will be presented 
where such decays occur with the desired strength. 

It is crucial that the non-standard decays of the muon are lepton number violating 



14| . Attempts to explain the LSND data in terms of lepton number conserving four- 
fermion interaction [JT5], [IToj will lead to inconsistencies with other experiments. One can 
make an SU(2)l transformation on such effective operators which would invariably gen- 
erate lepton number violating processes such as fi — > 3e decay, muonium-antimuonium 



transition, r — > \xee decay, etc, which are highly constrained by experiments |L4|, |T5j| . In 
contrast, the AL = 2 processes that we find will not admit an SU(2)l transformation to 
generate one of these highly forbidden processes, primarily because there is no charged 
current counterpart to the anomalous muon decays. 1 Since the new physics affects only 
^ decays, and not that of n ± , our interpretation of the LSND results would predict that 
no oscillation signal should be seen at the Fermilab mini-BOONE experiment. We have 
identified two AL = 2 effective operators which are invariant under the Standard Model 
symmetries that can induce the decay of Eq. (1) without generating any other effects. 
These effective operators can be realized from underlying renormalizable theories which 
are extensions of the Standard Model involving additional scalar multiplets. The masses 
of these scalar bosons are bounded by about 500 GeV, in order that the signal at LSND 
is significant. Since the scale of new physics, A, is relatively low, the particles at the scale 
A do not decouple entirely and have an impact on low energy observables. We highlight 
the processes that would succumb to the new physics most easily. They include small 
but observable deviations in e + e~ — > and e + e~ — > vvq cross sections, observable 

corrections to the muon g — 2, and lepton number violating Z° decays with branching 
ratios of order 10 -7 . 



2 Lepton number violating muon decay 

The decay fi + — > e + + V e + Vi with i = e fi or r, if allowed by the theory, can explain the 
LSND neutrino anomaly. Recall that the usual \i decay fi + — >• e + + z7 M + v e has no V e 
which the LSND detector registers. The new decay of /i + produces F e , which will then 
be detected. The branching ratio for the new decay should be about (1.5 — 3) x 10~ 3 , 
as can be inferred from the LSND analysis of the V^ — V e oscillation probability. 

We have found two effective operators invariant under the symmetries of the Standard 



Model ]T7| that can induce the desired AL = 2 muon decays. They are 

o 1 = ^(^e^x^cr 1 ^^) 



: We shall see shortly that the anomalous muon decay is mediated by neutral current. 
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e) L , etc denote the left-handed lepton doublets, while 



Here ^ 

eR, denote the right-handed singlets. $ is the Standard Model Higgs doublet with 
its hypercharge normalized to be +1/2. 

Note that both 0\ and O2 are non-renormalizable operators of dimension 9, and 
hence are suppressed by fifth power of A which characterizes the scale of new physics. 
As we shall see, these operators will arise from integrating out scalar fields which have 
lepton number violating interactions. Operator 0\ has a unique SU(2)l contraction, 
owing to Bose symmetry in the $ field. Although (^e^) in 0\ contracts to form an 
SU(2)l singlet, it will not induce a /7 L e^ mass term (see below). Similarly, no d = 5 
neutrino mass term will arise from these operators. Note also that there is a unique 
Lorentz contraction of fermionic fields in 0\, if we limit ourselves to obtaining these 
operators by integrating out scalar fields. Similar remarks apply to operator O2 as well. 
In this case, there is a minor variant, obtained by contracting (j!^^) rather than 
(jI R ty e &), which is somewhat different for the case when i 7^ e. 

When the vacuum expectation value ($°) = v is inserted in Eq. (2), it would generate 
the four-fermion operators 
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(3) 



These operators will lead to the decay \i 
L-violating /i + decay is given by 



Br(fj,' 



e V P V; 



e + + v e + Vi. The branching ratio for the 



(4) 
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Inserting the numerical value of v ~ 246 GeV, and demanding that the branching ratio 
for this decay is in the range Br(fi + — > e + V e Vi) = (0.0015 — 0.0025), we find that the 
scale A ~ (360 — 340) GeV. Since A is not very large, the particles that have masses of 
order A will not entirely decouple and will affect low energy observables. 

An important point about operators 0\ and O2 of Eq. (2) is that they induce only 
the terms given in Eq. (3) and nothing more. This is because of the AL = 2 nature 
of these operators and the fact that there is only one Higgs scalar in the theory. This 
proves that the effective operators of Eq. (2) generate the desired interactions that can 
explain the LSND events, and nothing else. Unlike in the case of AL = muon decays, 
there is no conflict with processes such as \i — > 3e, r — > /xee, etc. 



KARMEN experiment |]T8[ has set severe limits on lepton number violating decays of 
the muon. These limits however, depend sensitively on the assumed decay distribution. 
The anomalous muon decay arising from Eq. (2) in our models has a different distribu- 
tion compared to that of the usual (V-A) decay of the muon. We find that operators 
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0\ and O2 lead to the prediction p = for the Michel parameter for the AL = 2 decay. 
(The other decay asymmetry parameters for the L-violating p, decay are found to be: 
r\ — 0, £ = —3/4, and 5 = [0].) The limit from KARMEN for the branching ratio for 



the decay p + — > e+F e +Fj corresponding to the case where p = is Br < (0.0015 — 0.002) 
ipOfl . This limit is somewhat weaker than KARMEN's published limit corresponding to 
the case of p = 3/4 []2T|. Thus, the KARMEN results appear to be just about consistent 
with LSND observations, in our interpretation of the data. A joint analysis of the LSND 
and the KARMEN data in our framework will be desirable. 

All new phenomena in the model will arise from particles with masses of order A. 
To see these additional new effects, we have to generate operators of Eq. (2) from 
renormalizable Lagrangian densities. We now turn to this task. 



3 Gauge models for Lepton number violating muon 

decays 

We shall now present renormalizable gauge models that induce the effective operators 
of Eq. (2). These models are obtained by extending the scalar sector of the Standard 
Model. The new interaction Lagrangian for the first model which induces 0\ is given 
by 

+ M H^ X ® + \'A X && + H.C. (5) 

Here A(3, 1), x(3,0) and 1/2) are scalar fields (their SU{2)l x U(1)y quantum 
numbers are as indicated) which do not acquire any vacuum expectation value. 

Figure 1 will induce the effective operator 0\ from Eq. (5). Notice that the in- 
termediate scalar particles are electrically neutral, so the anomalous muon decay is a 
neutral current process. It is not possible to make an SU(2)l transformation in Fig. 1 
to convert some of the external fermions, without also transforming the external Higgs 
field to its charged component, which of course is unphysical. 

The strength of the effective four-fermion coupling arising from Fig. 1, in the ap- 
proximation of small scalar mixing, is 

n h^ e f ei X'M v 3 

Lr eS — ~ 2 9 2~T ' V 3 ) 

{m^ m z HO m z AO ) 

If all the mass parameters are equal and if h^ e , f e i and A' are equal to one, LSND data 
would require the scalar masses to be in the range (340—360) GeV. Actually, we can make 
this statement more precise. Let us denote the scalar mass-squared matrix involving 
the mixing of (H°,x°, A *) to be M 2 . Let K be the unitary matrix that diagonalizes 
M 2 : K^M 2 K = diag.(M 2 } Mf, Mf) with M 1 < M 2 < M 3 . Then by making use of the 
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Figure 1: Scalar exchange inducing lepton number violating decay /i + — > e + + z/ e + z/j 
through operator 

unitarity of K and the fact that (A4 2 )i3 = 0, we can write down the branching ratio to 
be 

Br{^ + ^ e + v e Vi) = ^——— 2 1 - —=) 1 - . (7) 



32M?G 2 F \ Ml) \ M. 



Noting that ji^-K^I < 1/2, and demanding that Br is in the range (0.0015 — 0.002), we 
obtain an upper limit on the lightest neutral scalar mass M\. M\ < (442 — 412) GeV. If 
|JCi2-K"3 2 | is less than 1/2, M\ will have to be even lighter. The heavier masses M 2 and 
M 3 cannot be very much larger, or else the mixing angle ji^-fQ^I wm t> e suppressed. We 
estimate M 2 3 to be not larger than M\ by about a factor of 2. While these limits hold 
for the neutral scalars, their charged partners should also have comparable masses. Any 
splitting between the masses of the neutral and charged scalars will be SU (2) ^-breaking, 
and is limited by the electroweak p parameter. 

in Eq. (5) can be ty e , or \I/ T . Consider the case when = \l/ e . In this case, the 
Lagrangian of Eq. (5) conserves two combinations of leptons numbers, viz., (L e + 3L M ) 
and L T . A Z% subgroup of electron number is also preserved. It is not required that 
these symmetries be broken, so they may be used to prevent rare processes at undesirable 
rates.Q 

If = in Eq. (5), then muon and tau lepton numbers are unbroken. Also 
unbroken is a Z 2 subgroup of electron number. If = \l/ r , then (L e + 2L M ) and 
(Lip + L T ) are unbroken, as is Z 2 subgroup of L e . These symmetries guarantee that 
potentially dangerous lepton number violating processes remain small. 

The scalar fields A, x an d H carry lepton number symmetries (or a discrete subgroup 
of these symmetries), so that they do not acquire VEVs. This can be ascertained by 
choosing the mass-squared of these scalars to be positive. The lepton number symmetries 

2 We have in mind a scenario where the solar and the atmospheric neutrino oscillations arise from 
neutrino masses and mixings induced by the seesaw mechanism. Lepton flavor violation that arises 
from the seesaw mechanism is too tiny to be observable in any processes other than neutrino oscillations 
themselves. 



will forbid possible tadpole contributions to their VEVs. Thus, these new interactions 
do not induce neutrino masses at all. The only source of neutrino mass in these models 
is through the seesaw mechanism. 

One may worry about the compatibility of the lepton number symmetries present in 
Eq. (5) and the neutrino oscillation data, which calls for the breaking of all such sym- 
metries. We shall present a concrete example to demonstrate the consistency. Consider 
Eq. (5) with = ty^. Let us impose L M and L T symmetries as well as a Z 2 sub- 
group of L e . Under these symmetries the scalar fields (if, x? A) have charges (1,0)_, 
(1,0)_ and (—1,0)_ respectively. The charged lepton mass matrix as well as the Dirac 
neutrino mass matrix will be diagonal due to these symmetries. Neutrino mixings can 
still arise from the superheavy Majorana mass matrix of the Vr fields. The charges of 
{veRi v iiR-, v tr) under these symmetries are {(0, 0)_; (1,0) + ; (0,1)+} respectively. If 
scalar fields with charges (0,0)+, (0, — 1)_, (—1,-1)+ and (0, — 2)+ are introduced and 
given large vaccum expectation values, all neutrino flavors will mix with one another. 
These superheavy scalars will have no couplings to the (H, \i A.) fields and thus will 
not spoil the symmetry, except through small neutrino mass terms. 

The exchange of the new scalars with masses of order 500 GeV can lead to new 
processes. We list below the most significant of these. 

1. The effective Michel parameter p in p decay is modified in our scenario: (p = 
0.7485). The deviation from 3/4 arises because of the rare p decay mode we 
have introduced here to explain the LSND data. Currently the uncertainty in 
the measurement of p is ±0.0026, but the TWIST experiment at TRIUMF |22| 



will have a sensitivity of 10~ 4 in p, which can test this scenario. There will also 
be a small shift in Gf extracted from p decay by (0.15 — 0.2)%, compared to 
the Standard Model value. Such a shift is currently consistent with radiative 
corrections in muon decay (parametrized by Ar), which has an uncertainty of 
about 0.2% arising from the ±5 GeV uncertainty in the top quark mass alone. 

2. The neutral component H° from the scalar H can mediate the process e + e~ —> 
p + p~ . The total cross section as well as the forward-backward asymmetry mea- 
sured away from the Z° pole give useful constraints. The most stringent one arises 

from LEP experiments run at 130 < y/s < 189 GeV. L3 collaboration has quoted 

o 



lower limits on contact interaction [23 at these energies. The exchange of H 



induces the following effective contact interaction: 

\h el 2 _ 

£ = -Tr^-{^Lln^L)e R ^e R ) . (8) 

This corresponds to the case of 7]l R = — 1 and all other rjij = in the notation of 
Ref. [23]. The limit on compositeness scale A„ > 1.9 TeV [Z3] implies 



m H o/\ VI > 379 GeV . (9) 
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The constraints from forward-backward asymmetry measurements, which have 
uncertainty of a few % will be satisfied once Eq. (9) is met. 

The exchange of A ++ partner of A will also contribute to the process e + e~ — > 
, The effective Lagrangian for this process is 



£ = 7rT—^e L )(l lL l%L) ■ (io) 



The constraint on A + is A+ > {3.8, 7.3, 3.9} TeV for i % = (e, fi, r) |2§. This 
would lead to the following limits: 

_^ > 7 5 8 GeV, I^f> 1456 GeV, > 778 GeV , (11) 

I Jee | | Je/i | | Jer | 

corresponding to ii being e, fi, or r. Since A ++ cannot be split much in mass 
from A , explaining the LSND result at the suggested rate would require f ei to be 
somewhat smaller than one and simultaneously one of the neutral scalar masses 
to be lower than about 300 GeV (see Eqs. (6) and (7)). 

3. The charged member H + from H scalar will mediate the process e + e~ — > z/z/7. This 
process has been studied at LEP as a way to count the number of light neutrino 
species. L3 collaboration quotes N v = 3.05 ± 0.11 ± 0.04 from single photon 



measurements carried out at 130 GeV < ^/s < 189 GeV |24j]. (Measurements 
very close to the Z° pole are less useful for our purpose.) Using the 3 sigma limit 
we obtain (following the procedure outlined in Ref. 



"Wl VI > 375 GeV ■ ( 12 ) 

We see that the predicted deviation from the Standard model is in the observable 
range, once the LSND data is explained. 

4. The neutral H° will contribute to the anomalous magnetic moment of the muon. 
The shift in a M compared to the Standard model prediction is given by 5a M ~ 
\h lie \ 2 (m fl /m H o) 2 /(2A7r 2 ) ~ 47 x 10" 10 | Vl 2 ( 100 GeV/m H o) 2 . This is in the exper- 
imentally interesting range for \hfj, e \ being order one and rriH° = (100 — 300) GeV, 
as needed by the LSND data. 

5. A + can mediate e + e _ — > vw-f for the case where = ty e , with a branching ratio 
very close to the current limits. In the case where \&, = there is new contribu- 
tion to ordinary muon decay from the exchange of A + . We obtain m&+ /|/ e ^| > 525 
GeV using 5(Ar) < 0.003. This limit, for the case where = ^ r At , while consistent, 
will force some of the other scalars to be lighter than 400 GeV. 

6. At a future linear collider running in the e~e~ mode, it is possible to produce the 
doubly charged scalar A ++ as an s-channel resonance in the process e~ + e~ — ► 
£~£~, which will provide a spectacular signal. 



7 



7. In the dimension 9 operator of Eq. (4), if we insert VEVs to the external Higgs 
fields, it becomes a d — 6 operator. Then we can attach a Z° boson line on 
the effective four-fermion operator. That would lead to rare Z° decays such as 
Z° — > e + [i^vv. A quick estimate gives the branching ratio for this decay to be of 
order (10~ 7 — 10 -8 ), which may be observable. 

8. The Z° boson can decay into a virtual HH pair, which can lead to Z° — > e + e~fi + fi~ 
signal. This branching ratio is smaller than 10~ 6 . 

3.1 Model for Operator 2 

The gauge model that induces operator C 2 in Eq. (2) can be obtained in analogous 
fashion. The Lagrangian of this model is taken to be 

+ MoH^x® + \'Ax&& + H.C. (13) 

The main difference from C\ is that the helicities of e and /i have been switched. The 
diagram of Fig. 2 will induce the decay ji + — > e + + v e + 77j. The phenomeno logical 
implications of the model are very similar to Model 1. The only significant difference 
is that the exchange of H + now cannot induce the process e + e~ — > uu^f, and so the 
constraint from that process would not apply. 




Figure 2: Scalar exchange inducing lepton number violating decay fi + —> e + + u e + z/j 
through operator 2 . 



4 Conclusions 

We have presented a class of models where the LSND neutrino anomaly can be explained 
in terms of the lepton number violating decay /i + — ► e + + u e + z/j. There are two effective 
operators that can lead to this decay. These operators do not cause any problem with 
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other rare processes such as \i — > 3e, unlike the case of AL = decays. We have presented 
gauge models where these operators are derived from renormalizable interactions. All 
models predict the existence of additional scalar fields with masses below about 500 
GeV. Observable deviations in the scattering processes e + e~ /U + /i~, e + e~ — ► i/i/j and 
lepton number violating Z° decays are expected with strengths not much smaller than 
the current experimental limits. 

Mini-BOONE experiment at Fermilab is expected to go on line in the near future. 
Our scenario predicts that mini-BOONE should see no signal. This is because the new 
interactions only affect ^ decays in our scheme, and not 7r =1= decays. Hence, a null 
result in Mini-Boone does not invalidate the LSND results. The proposed ORLaND 
experiment [^6| or experiments at a neutrino factory [TZTj using neutrinos from \x decays 
|l9fl can confirm or rule out the existence of such a lepton number violating decay 
mode; with the additional prediction that there be no L/E dependence. The TWIST 
experiment at TRIUM F |22|] can test for the shift in the effective Michel parameter that 
is predicted in our model. 
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